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An electric propulsion thrust stand capable of supporting testing of thrusters having a total mass of up to 1 25 
kg and producing thrust levels between 100 /zN to 1 N has been developed and tested. The design features 
a conventional hanging pendulum arm attached to a balance mechanism that converts horizontal deflections 
produced by the operating thruster into amplified vertical motion of a secondary arm. Tire level of amplification 
is changed through adjustment of the location of one of the pivot points linking the system. Response of the 
system depends on the relative magnitudes of the restoring moments applied by the displaced thruster mass 
and the twisting torsional pivots connecting the members of the balance mechanism. Displacement is measured 
using a non-contact, optical linear gap displacement transducer and balance oscillatory motion is attenuated 
using a passive, eddy-current damper. The thrust stand employs an automated leveling and thermal control 
system. Pools of liquid gallium are used to deliver power to the thruster without using solid wire conenctions, 
which can exert undesirable time- varying forces on the balance. These systems serve to eliminate sources of 
“zero-drift” that can occur as the stand thermally or mechanically shifts during the course of an experiment. 
An in-situ calibration rig allows for steady-state calibration before, dur ing and after thruster operation. Thrust 
measurements were carried out on a cylindrical Hall thruster that produces mN-level thrust. The measurements 
were very repeatable, producing results that compare favorably with previously published performance data, but 
with considerably smaller uncertainty. 

PACS numbers: 


I. INTRODUCTION 

Electric propulsion (EP) systems provide high specific 
impulse, but low thrust, relative to chemical propulsion 
systems. While chemical rocket thrust is generally mea- 
sured using load cellsfl], the low thrust levels associ- 
ated with EP devices lead to thrust stands which more 
closely resemble sensitive laboratory mass scales, where 
physical displacement of one of the mechanical mem- 
bers in the balance system is used to infer the applied 
force. These balances are typically configured as either a 
hanging, inverted or torsional pendulum with the thruster 
mounted at the end of the pendulum arm. 

The conventional hanging pendulum is the most sim- 
ple of the three configurations and is highly stable when 
subjected to external perturbations. However, high sen- 
sitivity can only be attained with a long pendulum arm, 
which may be impractical if the test facility (vacuum 
chamber) is small. Therefore, conventional hanging pen- 
dulum thrust stands have been primarily used to test high 
thrust-to-weight, T/W, EP devices such as electrother- 
mal arcjets and microwave electrothermal thrusters[2, 3], 
Also, the displacement sensor is typically mounted on a 
separate reference structure that can (undesirably) move, 
independent of the pendulum mounting structure, and 
produce erroneous results. 

Inverted pendulum thrust stands are less stable but 
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possess a greater sensitivity relative to the hanging pen- 
dulum configuration. They have been widely imple- 
mented to measure the performance of thrusters over a 
broad range of T jW, including resistojets, electrother- 
mal arcjets, MFD thrusters [4-6], Hall thrusters[7-9], 
and ion thrusters. A problem for the inverted pendulum 
configuration is that its stability is a strong function of the 
stiffness of the (supporting) flexures, which can change 
during the course of a test due to, for example, heating. 

Unlike both the conventional and inverted pendulum 
configurations, the restoring force in the torsional pen- 
dulum configuration (which rotates about an axis that 
is parallel to the gravity vector) can be made indepen- 
dent of the thruster mass. Consequently, torsional pen- 
dulum thrust stands provide the highest sensitivity, and 
have recently been implemented in micropropulsion per- 
formance evaluation [10-16]. They have also been shown 
to be effective in measuring the performance of pulsed 
thrusters[ 17-20]. The primary disadvantages of the tor- 
sional pendulum configuration are that the horizontal, 
asymmetric arrangement of its members can be difficult 
to configure within a vacuum chamber possessing lim- 
ited space and its sensitivity does not readily allow for 
testing of higher power thrusters. 

We describe in this work a new thrust stand design: the 
Variable Amplitude Hanging Pendulum with Extended 
Range (VAHPER)[21 ]. This thrust stand is a variant 
of the conventional pendulum, but it incorporates many 
novel features designed to enhance the quality of thrust 
data and the range of devices that can be tested. Specif- 
ically, the VAHPER thrust stand was designed to en- 
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balance mechanism 



FIG. 1; Schematic illustration of the assembled VAHPER 
thrust stand. 


described in this paper. Data demonstrating the sensitiv- 
ity and thrust measurement capabilities of the VAHPER 
thrust stand are presented in section IV 


II. VAHPER SYSTEM OVERVIEW 

A schematic illustration of the VAHPER thrust stand is 
shown in Fig. 1. In die rest of this section, we describe in 
detail the various subsystems and die enhancement each 
component adds to the fidelity of the thrust measurement. 


A. Balance Mechanism 

The system of linkages comprising die VAHPER 
thrust stand balance mechanism are represented as a sim- 
plified schematic in Fig. 2. All of die pivots (points 
B,C,D,E) utilize frictionless torsional flexures manufac- 
tured by the C-Flex Bearing Corp. Member ABC piv- 
ots rigidly about point B while member DEF pivots 
about point E. Rotation of the two members are coupled 
through linkage CD. Unless otherwise noted, displace- 
ments are represented by the vector .dy, member lengths 
by Ijj and angles by flyt where i, j, and k are the end- 
points. 

Application of thrust T causes members ABC and 
DEF to rotate through angles a% and a^, respectively, 
to the new position indicated by the primed (') points. 
In this configuration, horizontal deflection d A A' is con- 


able thrust measurement of electric thrusters weighing 
between 1-125 kg and producing thrust levels from 1 mN 
to over 1 N. This is accomplished through a mechanical 
linkage system that translates horizontal displacement 
into amplified vertical displacement. The system pro- 
vides for adjustable sensitivity though the adjustment of 
a single member. Displacement monitoring is performed 
on tire vertically displaced member, which eliminates the 
need for a separate reference structure. 

Pendulum-based thrust balances are susceptible to 
‘zero-drift’, which means that when the thruster is turned 
off, the pendulum arm does not return to its initial posi- 
tion. Several components and subsystems have been in- 
cluded in the VAHPER thrust stand design to mitigate the 
causes of zero-drift. These include low friction flexural 
pivot mounts, a non-contact ratiometric position sensor, 
an on-board level control system, a thermal management 
system and an innovative method to feed power onto the 
thrust stand arm using liquid gallium pools. 

In the next section we provide details regarding the de- 
sign of the thrust stand. Section III contains information 
about the thruster and facility used for the experiments 



FIG. 2: Schematic illustration of the balance mechanism com- 
ponents in their undeflected and deflected positions (deflected 
beam positions shown in grey). 
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verted to (amplified) vertical deflection dpF'. which is 
measured by the displacement transducer at the end 
of beam DEF. The VAHPER thrust stand pivots are 
mounted to the same structure as the displacement trans- 
ducer, thus eliminating the need for a separate (and prob- 
lematic) reference structure. The trigonometric exercise 
of relating <1 aa' to dpr' is performed in Appendix 1 and 
results in the equation set 



Oi = <*1 + 0BC'D' + ^C'D'E “ ^BCD — ^CDE, (1) 


where all the lengths and angles are known in terms 
of the initial geometric configuration and are defined in 
Appendix 1. The coordinates of the endpoints in the 
VAHPER thrust balance are given in Table I. Given 
this geometry, the amount of displacement amplification 
(|dFF'| / Maa'I) realized in the system is plotted in in 
Fig. 3 as a function of the variable length Zde- 

To find the amount of actual deflection |cIff' I realized 
for a given level of thrust, we must perform a free body 
analysis of the various members, taking into account the 


TABLE I: Coordinates of the undeflected points comprising the 
VAHPER thrust stand as shown in Fig. 2. 



x [cm] 

V [cm] 

A 

0.0 

-121.9 

B 

0.0 

0.0 

C 

31.4 

-3.2 

D 

0.5 to 3.0 

4.2 

E 

30.4 

4.2 

F 

-12.9 

4.2 



FIG. 3: Displacement amplification in the VAHPER thrust bal- 
ance as a funcdon of the variable length fo B . 


TABLE H: Spring constants for the torsional pivots in the VAH- 
PER thrust stand. 



Spring Consumt [cm-N/rad] 

kB 

153.4 

he 

4.21 

/CD 

4.21 

kB 

8.42 


restoring moments exerted by the torsional flextures and 
the mass of the thruster (located at point A). This analy- 
sis is found in Appendix 2 and results in an expression 
relating the thrust |T| to the angle a\'. 


T = (M b + Me ) /l ab + m g sin (qj. ) 


(Mg + Mp>) Zbc cos | 

GHbcmv) 

/aB /de COS - 

- #C'D'E j 


( 2 ) 


where the torsional pivot reaction moments M are de- 
fined in Eqs. (A12) with spiing constants given in Ta- 
ble II. After the variable length Ip E is fixed, every term 
on the right hand side of tire equation is a function of 
only cti. The problem is most easily solved by iteratively 
varying c*i until Eq. (2) equals a given thrust level. Equa- 
tions (I) can then be used to find the member deflections 
in the system. In Fig. 4A we have plotted the predicted 
deflection at point F as a function of l de for a thruster 
having a mass of 5 kg and producing a fixed thrust of 10 
mN. The same type of plot is shown in Fig. 4B for the 
case where a 40 kg thruster produces a constant 20 mN 
of thrust. 

From Fig. 4, we see that even though the displacement 
amplification increases rapidly with decreasing Zde. the 
actual deflection produced at point F by a given thrust 
level does not necessarily increase over the same range 
of Zde but, rather, an optimal value exists for a given 
thruster and choice of springs. The reason for this is 
clearly illustrated in Fig. 4A, where moments exerted by 
the torsional pivots in response to the motion induced 
by the thruster begin to dominate the thrust stand’s re- 
sponse as we decrease /de below 1 .5 cm. Even in Fig. 
4B where the thruster is significantly heavier, we observe 
the deflection curve beginning to peak as l de approaches 
5 mm, indicating that the restoring forces in the system 
arising from the reactions of the torsional pivots are be- 
ginning approach the same level as the restoring force 
applied by the thruster’s mass. While in this paper we 
only focus on cases that are represen tative of the existing 
VAHPER thrust stand hardware, a parametric analysis 
using Eqs. (1) and (2) could be performed to determine 
the global sensitivity of the curves in Fig. 4 to the various 
member lengths and spring constant values. 
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FIG. 4: Predicted deflection, j, as a function of the vari- 
able length Ibe for A) a 5 kg thruster producing 10 mN of thrust 
and B) a 40 kg thruster producing 20 raN of thrust (spring con- 
stant values from Table II are assumed in both cases). 



FIG. 5: Minimum and maximum measurable thrust levels on 
the VAHPER stand as a function of thruster mass and the vari- 
able length Ide- The minimum and maximum levels corre- 
spond to deflections at point F of 13 /rm and 8.5 mm, respec- 
tivley (spring constant values from Table II are assumed in both 
cases). 


B, Displacement Sensing 

Measurement of the distance |cIff'| is achieved us- 
ing a BEI Precision Systems Linear Gap Displacement 
Transducer (LGDT), which is a non-intrusive (not in 
physical contact with member DEF) optical transducer 
that is relatively immune to external electric, radio fre- 
quency (RF), trad magnetic interference. The LGDT il- 
luminates the target surface with an IR LED through a 
concentric-ring fiber-optic probe. Light reflected from 
the target is analyzed using a ratiometric technique to 
derive the distance to the target surface. The analog volt- 
age output of the LGDT varies linearly with the distance 
from the target to the tip of the probe. 

The LGDT used on the VAHPER thrust stand has a 
resolution of 1.3 pm, a sensitivity of 0.75 mV/pm, and a 
1 % linear range of 8.5 mm. The sensor output is trans- 
mitted through a single, shielded twisted pair to a 12-bit 
A/D converter (NI PXI-1601). This conversion brings 
the overall resolution of the system to about 13 /tm. Tak- 
ing this as the minimum measurable distance |cIff'| and 
8.5 mm as the maximum, and using Eq. (2), we can solve 
for the minimum and maximum resolvable thrust levels 
as a function of thruster mass and the variable length 
Idb- These bounds are plotted in Fig. 5 for the VAH- 
PER stand. We observe that for Jow thruster masses 
where the reactions of the torsional pivots dominate the 
thrust stand response, the minimum resolvable thrust is 
measured when the stand in in the configuration yielding 
the smallest displacement amplification (greatest value 
of (de)- in addition, the greatest measurable thrust level 
is obtainable in the configuration yielding the largest dis- 
placement amplification (smallest value of (de)- How- 
ever, as thruster mass is increased, this counterintuitive 
trend reverses and the configuration yielding the greatest 
displacement amplification allows for measurement of 
the minimum resolvable thrust levels while the configu- 
ration giving minimal displacement amplification should 
be employed to measure the maximum resolvable thrust. 


C. Damping 

Attenuation of mechanical oscillations in the motion 
of the thrust balance is achieved using a passive, eddy- 
current damper. The damper is located at the end of the 
balance secondary beam, near the displacement sensor 
(see Fig. 1). The damper assembly is pictured in Fig. 6 
and is comprised of four 3/8”-diameter samarium-cobalt 
permanent magnets embedded in a soft iron yoke. Op- 
posing magnets are separated by a 1/8” slot, providing 
a magnetic field strength in the gap of —0.2 T. A 1/16”- 
thick copper plate, attached to the end of the secondary 
beam, extends into the slot between the permanent mag- 
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FIG. 6: Eddy-current oscillation damper illustration (note: iron 
yoke rendered transparent to reveal inner detail). 


nets. Vertical motion of the secondary beam causes the 
copper plate to move relative to the applied magnetic 
field, inducing eddy-currents in the copper that dissipate 
the mechanical energy and provide near critical damp- 
ing of thrust stand oscillations. The amount of dissipa- 
tion can be adjusted either by adjusting the magnetic field 
strength or by changing the size or material of the mov- 
ing plate to alter the exposed area or conductivity, respec- 
tively. 


D. Level Control 

Zero drift of the displacement sensor signal is a ubiq- 
uitous problem in thrust measurement systems, deriving 
primarily from unsteady thermal loads and vacuum tank 
distortion during the course of a test. We have included a 
level control system in the VAHPER thrust stand design, 
to both sense and counteract deflections of the structure 
that supports the balance mechanism, with the aim of 
minimizing the zero drift dui'ing tests. 

Level sensing is performed using a pair of high ac- 
curacy Jewell Instruments inclinometers with 0.001° of 
accuracy over their full scale range of ±3° relative the 
the Earth’s gravity vector. The sensors are mounted to 
the baseplate (see Fig. 1) perpendicular to each other to 
allow for two-axis (X, Y) detection. Any deviation from 
0° (or any other user defined angle) causes the control 
system to adjust the baseplate position. The adjustments 
are performed by two high resolution (1 pm/step) four- 
phase stepper motors, which change the position of the 
baseplate until the inclinometers indicate that the base- 
plate is level. 


E. Thermal Management 

A closed-loop thermal management system was im- 
plemented, with the goal of maintaining a constant bal- 
ance temperature during tests, in order to minimize dis- 
placement zero drift. The balance mechanism is mounted 
on a thick (1.5”), water-cooled aluminum plate (see Fig. 
1) to provide high thermal inertia. The water tempera- 
ture is maintained at a constant temperature by a Thermo 
Neslab M150 programmable water chiller, which is ca- 
pable of removing up to 5 kW of heat power. Tempera- 
ture is measured on the thrust stand at various locations 
by thermocouples. 


F. Thruster Power Feeds 

In many thrust measurements, the wires leading to the 
thruster can greatly affect the level of zero drift expe- 
rienced over the course of a test. The wires experience 
Joule heating, causing them to deform during thruster op- 
eration (especially at high current levels like those used 
to pre-heat an emitter cathode), which can apply time- 
varying forces and torques to the thrust stand arm. In the 
VAHPER thrust stand design, this source of drift is elim- 
inated by conducting current onto the thrust stand arm 
through pools of liquid gallium. Since the wires on the 
arm are only in contact with the liquid metal, they are 
unable to apply loads to the thrust stand arm. 


G. In-Situ Calibration 

Displacement (thrust) calibration was accomplished 
by applying a series of known loads normal to the pendu- 
lum arm, and recording the resulting displacement mea- 
sured at point F. A set of five known masses is attached at 
regular - intervals along a Kevlar - string (see Fig. 1). The 
string is wound around a spool and attached to a stepper 
motor, which allows for remote actuation of the calibra- 
tion procedure. The bottom end of the string is draped 
over a low-friction pulley and attached to the hanging 
pendulum arm. As the weights are lowered, the string- 
mass-pulley system transfers the gravitational force on 
the weights to the thrust stand, normal to the pendulum 
arm. The calibration procedure is automated so that at 
any time (before, during or after a thrust measurement), 
the calibration weights can be lowered and retracted to 
provide continuous checks on the displacement calibra- 
tion. 
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FIG. 7: VAHPER thrust stand control system architecture and 
data flow. 


H. Data Acquisition and Control 

The control system and data flow are shown schemat- 
ically in Fig. 7. The system architecture is based around 
a National Instruments PXI-8176 embedded computer, 
which is mounted in a chassis along with several other 
modules that allow for communication with and control 
of various subsystems. Temperature and displacement 
data (in the form of a voltage signal) are brought into 
the system through a TX-68 terminal board. The step- 
per motors are controlled and by the M1D-7604 mod- 
ule, which also receives feedback data from the incli- 
nometers. Two-way communication with the program- 
mable water chiller is performed through an RS-232 con- 
nection, By communicating with other PCs thought an 
Ethernet connection, a common time stamp can be es- 
tablished. This ability allows for all data to be tagged 
with the same time stamp, regardless of which computer 
is acquiring data. This is critical in the post analysis 
phase when thruster data (current, voltage, mass flow 
rate), which is controlled by a separate system, is cross- 
referenced and correlated with thrust stand deflection. 


III. EXPERIMENTAL SETUP 

We conducted experiments using a low-power 
(0(100) W) cylindrical Hall thruster to test the abil- 
ity to resolve thrust levels on the lower end of the pre- 
dicted thrust stand range. For the purpose of validation, 
a thruster was chosen for which thrust stand measure- 
ments from another facility already existed. Based on 
those measurements, we expected thrust levels between 
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Region Neutralizer 

FIG. 8: A) Schematic illustration and B) photograph of the 
PPPL cylindrical Hall thruster. 


1-10 mN. 


A. Vacuum Facility 

The thrust stand is mounted inside a 9-ft diameter, 25- 
ft long stainless steel vacuum chamber. The vacuum 
level inside the chamber is maintained by two 2400 I/s 
turbopumps and two 9500 1/s cryopumps. The base pres- 
sure was 5.7 x 10~ 7 Torr and the background pressure 
of xenon during thruster operation at ~5 seem total flow 
rate was ~9 x 10~ 6 Tom 


B. Cylindrical Hall Thruster 

The Princeton Plasma Physics Laboratory (PPPL) 3 
cm cylindrical Hall thruster (CHT) was used for these 
experiments. A The CHT is illustrated schematically 
in Fig. 8A, while a photograph is shown in Fig. 8B. 
The Boron-Nitride thruster channel is a composite of a 
shorter, annular region and a longer, cylindrical region. 
Gas is injected through holes in the annular anode. The 
magnetic circuit consists of a magnetic core and two 
electromagnet coils that can be operated using indepen- 
dently controlled power supplies. More details on the 
CHT can be found in Refs. [22-26] and the references 
therein. 

The working propellant for these experiments was re- 
search grade xenon gas. The cathode and anode flowrates 
were independently controlled using two MKS 0-10 
seem flow controllers (calibrated on Xe). A commercial 
HeatWave plasma source was used as the cathode neu- 
tralizer. In all our experiments, the cathode flowrate was 
1 seem. 
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IV. EXPERIMENTAL RESULTS 

A. Thrust Stand Calibration 

Displacement (thrust) calibration of the VAHPER 
thrust stand is accomplished using an in-situ calibration 
rig discussed in the previous section. A displacement 
waveform, as shown in Fig. 9A, is obtained by applying 
a series of known loads normal to the pendulum arm. The 
stability in the displacement zero is evident as the stand 
returns to its initial position roughly 14 minutes after we 
initiate the calibration procedure. 

Calibration measurements and the resulting fit to the 
data set are plotted in Fig. 9B. An automated data reduc- 
tion routine is used to analyze the displacement data and 
generate the curve fit coefficients which can be used to 
compute thrust. This analysis is performed under the as- 
sumption that the relationship between the applied force 
and the measured displacement is linear. The error prop- 
agation algorithm, which perform a linear curve fit of the 
data taking into account the error bars on the individ- 
ual force and displacement measurements, is discussed 
in Ref. [27], Based on the curve fit, we see that the ulti- 
mate resolution without the thruster running is ~50 //N. 



LGDT Output [V] 


B. Raw Thrust Data 

The output of the LGDT, power supplies and mass 
flow controllers are recorded throughout the duration of 
the test. Raw data output by the LGDT is displayed in 
Fig. 10 for two separate CHT discharge voltages. We 
immediately see from this figure that there is very little 
zero-drift in the thrust balance position over this time pe- 
riod. Typically, the largest zero-drifts in past experiments 
have occurred between the time when the cathode heater 
is turned on and the discharge is ignited. This is due 
to the large current level needed to heat the cathode to 
thermionically emitting temperatures (23 A for our Heat- 
Wave cathode). We see in the VAHPER thrust stand data 
that there is no significant drift associated with cathode 
heater operation, thus validating our liquid metal power 
conduction scheme. 

Prior to thrust measurements, the thruster is operated 
for approximately 30 minutes to condition the insulating 
surfaces in the discharge chamber. Thruster operation 
begins when the emitter cathode becomes hot enough to 
sustain a thermionic discharge. At this point, the dis- 
charge voltage is applied between the cathode and the an- 
ode. After the main discharge is activated, the magnetic 
field coils are energized and adjusted to the desired cur- 
rent levels. As expected, we observe fluctuations in the 
raw data trace as the field is adjusted to the final value 
since the actual thrust varies as the L-field is adjusted. 


FIG. 9: A) Time history of VAHPER calibration data showing 
LGDT measured output as the calibration weights were added 
and removed. B) Applied calibration force plotted as a function 
of LGDT response. The linear curve is a fit to the displayed 
data. Error bars on the datapoints and the fit coefficients repre- 
sent a 68% ( 1 cr) confidence interval. 


Once the field is set, no further adjustments are made, al- 
lowing for a 1 .5-2 minute steady-state period. After this 
time period, the discharge is extinguished and the field 
coils are deenergized in preparation for the next test. We 
allow enough time between thruster operating points to 
re-establish the LGDT output baseline (zero level). 


C. Thruster Performance Measurements 

The thrust (T) is found by first finding the LGDT out- 
put voltage difference between the steady-state portion of 
thruster operation and the zero level output immediately 
following thruster shutdown (see Fig. 10). These data 
are then converted to thrust using the calibration curve- 
fit constants found in Fig. 9B. Specific impulse ( I sp ) and 
anode efficiency ( 77 ) are defined as 

- JL - 7,2 

sp mgo ’ * 2 fnP ' 
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FIG. 10: Raw LGDT output data as a function of elapsed time 
for two thruster operating points. 


where P is the discharge power and m is the mass flow 
rate through the anode. The uncertainty levels on the per- 
formance data represent a 95% (2 cr) confidence interval. 

Thrust measurements corresponding to an anode mass 
flow rate of 3.9 seem are presented in Fig. 11 A, with 
the corresponding I sp values found in Fig. 1 IB and the 
anode efficiencies plotted in Fig. 1 1C. Our data span a 
range between 130 to 185 W in discharge power, with 
resulting thrust levels between 4.8 and 6.2 mN, 7 sp lev- 
els between 1275 and 1650 s, and anode efficiencies from 
23 to 27%. Also plotted in this figure are measurements 
from Ref. [28] performed at Princeton University on a 
swinging-gate style thrust stand using the same thruster. 
The agreement between this data set and our data is very 
good, and the error bars on our data are ~20% of those 
associated with the measurements performed at Prince- 
ton. 


V. CONCLUSIONS 

We have developed a hanging pendulum thmst stand, 
which employs an adjustable mechanical linkage system 
to convert the horizontal deflection of the pendulum arm 
to amplified vertical deflection of a secondary beam. Fur- 
thermore, the geometry allows for mounting of the sec- 
ondary beam and the deflection measurement system to 
the same baseplate, eliminating the need to perform mea- 
surements of deflection at the end of the pendulum using 
a separate reference structure. Deflection amplification is 
a function of geometry, and depends only on the variable 
length Ide- The system’s response, on the other hand, 
depends on the relative magntitudes of the moments ap- 
plied to the system by the displaced thruster mass and the 
deflected torsional pivots. The thrust stand is capable of 




Discharge Power [W] 

FIG. 1 1 : Performance measurements for a 3 cm cylindrical Hall 
thruster: A) Thrust, B) I sp , C) Anode Efficiency as a function 
of discharge power. The error bars represent a 95% (2cr) con- 
fidence interval. No error bars are displayed on the data from 
Ref. [28]. 


measuring thrust levels ranging from 100 /iN to 1 N, and 
can support thrusters of up to 125 kg total mass. 

Several supporting systems have been implemented 
to enhance the measurement accuracy: high-sensitivity, 
non-contact displacement sensing, an automated leveling 
system, an oscillation damping device, thermal control, 
and a novel means to deliver power to the thruster. These 
subsystems serve to eliminate the sources of zero-drift. 
In addition, an in-situ calibration rig enables displace- 
ment calibration at any time during a test. 

Experiments using a well-characterized thruster vali- 
dated the capability of the thrust stand to measure thrust, 
especially at low thrust levels. Our data were in good 
agreement with a data set obtained using a different, 
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well-characterized thrust stand, the error bars on our 
data set were only 20% of those found in the literature, 
demonstrating the effectiveness of the systems listed in 
the previous paragraph. 
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APPENDIX: THEORETICAL ANALYSIS OF BALANCE 
MECHANISM 

We proceed with a derivation of Eqs. 1. and 2 from first 
principles. Scalar quantities of the displacement vectors 
dy are represented by dy. All endpoints and vectors are 
the same as those defined in Fig. 2. 


1. Balance Mechansim Deflection Analysis 

To find the displacement amplification, we must relate 
the displacements cipp< and (Iaa'- From Figs. 12A and 
12B we see that isosceles triangles AA'B aud EFF' are 
used to write the displacements in terms of the angles a% 
and 02 : 


ai 

= 2 sin -1 ^ 

d AA' \ 

,2Zab ; ’ 

(Al) 

d-FF' 

= 2 Ief sin 

(?)• 

(A2) 


Writing an equation relating ripp- and c?aa' now comes 
down to an exercise in writing o 2 in terms of ay. 

Once the length (de is chosen, the initial separation 
distances and angles within the linkage mechanism are 
fully known. The bulk of our analysis will focus on the 
quadrilateral BCDE shown in its undefiected position in 
Fig. 12C. When the mechanism is deflected as shown in 
Fig. 12D, we can write (see Fig. 2) 

$EBC' ~ #ebc + a l 1 (A3) 

#D'EB = 0DEB - «2- (A4) 


The length Ice is written using the law of cosines as: 


O 1 ^9 

Ice — {Iqc + ^be ~ 2 Ibc (be cos (#ebc + «i)) 

(A5) 

With all the lengths in the triangles BC'E and C'D'E 
known, we again use the law of cosines to write: 


6>i 


BC'D' 


COS 


$C'D'] 


+ cos" 


COS 


I 2 4- I 2 
BC ~L t C'E 


' (be 


2 (bc Ice 
2 

CD 


I'cE+lr, 


/2 

f DE 


2 Ice (cd 

(de + (cd ~ (q'e 
2 (de (cd 


(A6) 

(A7) 


The geometry is such that the angles comprising 
quadrilaterals BCDE and BC'D'E must both sum to 27r 
radians. If we equate the angles from the two quadrilat- 
erals, substitute for the angles in BC'D'E using Eqs. (A3) 
and (A4), and rearrange the terms, we obtain 


C*2 = Ol + 0BC'D' + BC'D'E “ 0BCD ~ #CDE, (A8) 


Given an initial geometry and selecting a value for a i, all 
the angles on the right-hand side of Eq. (A8) are known. 
Equations (Al), (A2) and (A8) allow us to relate d A A' to 
dppi and are exactly those given as Eqs. (1). 


2. Linkage Response to Applied Thrust 

An operating thruster applies a force T to the thrust 
stand, deflecting the linkages to a new position denoted 
by primes (') in Fig. 2. This motion is opposed by the 
torque arising from the weight of tire thruster, mg, and 



FIG. 12: Schematic illustrations (not to scale) of various mem- 
bers in the VAHPER thrust balance. The endpoints, vectors and 
lengths are the same as those defined in Fig. 2. 
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FIG. 13: Free body diagrams (not to scale) showing the relevant 
forces and torsional spring moments applied to the (A) upper 
and (B) lower linkages in the VAHPER thrust balance when 
the system is deflected. The deflected positions are defined in 
Fig. 2. 

by the moments produced as the torsional springs flex 
from their neutral positions. To find the actual deflection 
4 aa' we must sum the moments applied to each linkage. 

Free body diagrams showing the moments applied by 
deflected torsional springs and external forces applying 


additional moments to both the upper and lower linkages 
are shown in Fig. 13. In the static case, summing the 
moments for linkage A'BC' about point B yields 

A'/b + Me — T Iab T 7ii (j l \ b sin (c*i) 

+ ^C'D' Ibc cos - 0bc'D'^ = o, (A9) 

while doing the same for linkage D'EF' gives 

Me + Md' — Pc'D' Ide cos — 0c'd'e) — 0. (A10) 

Solving Eq. (A10) for Few, substituting that result 
into Eq. (A9) and rearranging yields: 

T = (Mb + Me) /Iab + rri g sin (cm) 

(Me + Md') Ibc cos — Pbc'D') 

4 _ v/ / i-Ain 

‘ / 7T 

Iab Ide cos \ 2 ~ ^ C ' D ' E 

which is the same equation given in Eq. (2). The mo- 
ments in Eq. (Al 1) are 

Mb = o . 1 ha , Me = (0bc'D' - #bcd) ^Ovait) 

M e = 0.2 ks, Me = (#c'D'E - #cde) &d> 

where the k values represent the different torsional spring 
constants at each position. From the deflection analysis 
of the previous section, all the angles in Eqs. (All) and 
(A12) can be computed given an initial deflection cIaa' 
(or more simply by choosing a value of ai). 
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